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We investigate optical transitions of non-interacting electron systems consisting of two symmetric 
energy bands touching each other at the Fermi energy (e.g. graphene at half filling). Optical 
conductivity is obtained using Kubo formula at zero temperature. We show that for particles 
whose pseudospin direction is determined solely by the direction of their momentum, the optical 
conductivity has power law frequency dependence with the exponent (d — 2)/z where d is the 
dimension of the system and z is the dynamical exponent. According to our result two-dimensional 
systems with the above pseudospin characteristics always exhibit frequency-independent optical 
conductivity. 

PACS numbers: 78.67.Wj, 78.67.Pt, 03.65.Aa, 64.60. F- 
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I. INTRODUCTION 

Graphene is a two-dimensional material consisting of 
carbon atoms arranged in a honeycomb lattice. Since 
its discover;^, it attracted lots of attention due to its 
peculiar electronic and optical properties^. Graphene 
is a promising material in industrial applications and 
is also very important from theoretical point of view 
because its quasiparticles can be described as massless 
Dirac fermions. These particles possess pseudospin de- 
gree of freedom originating from the two atoms of the unit 
cell. One of the most fascinating properties of monolayer 
graphene is that the optical conductivity is universal and 
independent of the frequency^- 7 - . This feature was argued 
to be a consequence of the linear spectrum and the struc- 
ture of the Dirac cones^. 

Electronic and optical properties of multilayer 
graphene have also been investigated experimentally 
and theoreticallj*2r— . The band structure of multilayer 
graphene is very sensitive to its stacking sequence^. In 
the case of periodic ABC stacking, the low energy be- 
haviour of quasiparticles show chiral nature and can be 
characterized by the pseudospin winding number which 
is equal to the number of layers. Chiral nature of quasi- 
particles has been proven to play an important role in 
various phenomena, such as quantum Hall effect^, THz 
radiation measurements^ and angle-resolved photoemis- 
sion spectroscopy measurements^. 

In periodic ABC stacked (sometimes called chiral or 
orthorhombic) multilayer graphene the Hamiltonian can 
be written as 



conductivity of ABC stacked multilayer graphene is also 
found to be universal and independent of frequency^. 

In the present paper the optical conductivity of two- 
band systems is studied within Kubo formalism. We re- 
view earlier calculations^ and obtain the real part of the 
optical conductivity expressed in terms of the pseudospin 
of the particles. Analytical results are obtained for two- 
band systems in which the pseudospin of the quasiparti- 
cles is determined only by the direction of their momen- 
tum (as well as in the case of chiral multilayer graphene). 
The optical conductivity of general chiral systems defined 
in this way exhibits power-law frequency dependence. 
This behaviour is a consequence of the fact that there 
is no energy scale in these systems, therefore, scaling ar- 
guments can be applied. 

We show that the universal behaviour of the optical 
conductivity in chiral multilayer graphene can be ex- 
plained by the two-dimensionality and the invariance un- 
der dilatations. 



II. OPTICAL CONDUCTIVITY IN SOLIDS 

Physical quantities measured in optical experiments 
can usually be expressed by means of the optical conduc- 
tivity. Using the Kubo formula the optical conductivity 
can be written as 
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where k and 7 are the magnitude and the angle of the vec- 
tor k which is the wavenumber measured from the corners 
of the Brillouin zone. In the formula, M is the number of 
layers and v = ±1 is the valley index describing the two 
incquivalent corners of the Brillouin zone. The optical 



where n(r) is the density of electrons, j a (r, t) is the cur- 
rent density operator and the expectation value shall be 
evaluated in the equilibrium state. The first term in (JTTJ) 
is the diamagnctic term. Henceforth, we neglect the no- 
tation lim i 5_ > o+ • 
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The one-particle eigcnstates of a solid are the Bloch 
states which obey the i/(r)* /k (r) = Ei(k)^i k (r) 
Schrodingcr equation where H(r) = —h 2 A/(2m) + U(r) 
with lattice periodic potential. Because of Bloch's the- 
orem the eigenfunctions can be rewritten as 'l'ik(r) = 
e lkr ttik(r) / V?2 where uzk(r) is a periodic function of 
space variables and f2 is the volume of the sample. If 
the Coulomb interaction between electrons is neglected 
then one-particle eigenstates are filled following the Fermi 
distribution function f(E). 

The q = Fourier component of the frequency depen- 
dent conductivity 
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where k and k' run over the whole Brillouin zone. The 
matrix elements Vu<kk' = (^ikWl^i'k') of the velocity 
operator may be rewritten as 



V«/ kk / = S kk , [iuw (k)A u ,{k) + vi(k)5 u >} 
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with huiw (k) = Ei (k) — Ev (k) . In the formula, hvi (k) = 
dEi I dk is the quasiparticlc velocity in the Ith band and 
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is the self-adjoint Berry connection matrix. Note that the 
first term in Q describes transitions between two differ- 
ent bands and equals to zero if / = I'. The second term 
is non-zero only if I = I' and is the quasiparticlc velocity 
as expected. However, in ([3]) only terms corresponding 
to two different bands survive for non-zero frequencies 
because of the Fermi functions. Therefore, 
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In the followings, we focus on the real part of the op- 
tical conductivity. The imaginary part may be calcu- 
lated using Kramers-Kronig relations. The diamagnetic 
conductivity is purely imaginary so it does not appear 
in the real part at non-zero frequencies. The 5 — > + 
limit of the second term in ((5]) can be divided into a 
principal value part and a Dirac-delta part. Since the 
Hamiltonian does not depend on the spin, the system is 
time reversal invariant. This leads to Ei(k) = Ei(—k) 
and A;;/(k) = Aj/;(— k) and, hence, the real part of the 
conductivity equals to the Dirac-delta part. We note 
that this statement is true even if the Hamiltonian is 



spin-dependent, for example in the presence of magnetic 
field. In this case the statement is a consequence of 
Ef(k) = Ef s (-k) and A s w (k) = A~ s (-k) with s de- 
noting the different spin directions, but magnetic effects 
are out of scope of the present paper. Therefore, 
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can be obtained for non-zero frequencies. Note that 
intraband processes do not contribute for non-zero 
frequencies^. The optical conductivity is now expressed 
in terms of the Berry connection matrix. It is well known 
that the Berry connection is essential in determining the 
dc Hall conductivity^. However, Eq. Q shows that 
this topological quantity plays an important role in the 
optical conductivity as well. 

Using Eq. ([7]) one can calculate the optical conductiv- 
ity of solids generally. To do so, first the energy spectrum 
and the Berry connection matrix have to be computed in 
band structure calculations. In tight binding approxi- 
mation the Hilbert space is restricted to several atomic 
orbits per unit cell and the eigenfunctions are given as 
e tkr |k, i)/VO where |k, I) is a finite-dimensional vector 
representing the eigenstate on this subspace. In order 
to calculate the Berry connection matrix elements one 
may use Aw(k) = i(k, l\ V k |k, I') instead of (|SJ). This 
expression depends on the choice of the basis in the 
tight binding subspace but we claim that in most cases 
i(k, Z|V k |k, I') deviates from ([5]) negligibly in the "natu- 
ral" tight binding basis. For details, see Appendix. 



III. TWO BAND SYSTEM 

For frequencies lower than the bandwidth, only optical 
transitions between the two bands closest to the Fermi 
level are allowed. We focus on systems with two bands 
where the lower band is fully occupied and the upper 
band is empty and we consider zero temperature. 

The Hamiltonian of two band models can be repre- 
sented on this subspace with a self-adjoint 2x2 matrix. 
Assuming that it is independent of the spin the Hamilton 
operator can generally be written as 

H(k) = e {k)I + e(k)n(k)& 

where £o(k) and e(k) > arc real functions. The three- 
dimensional vector a is the pseudospin operator built up 
from the Pauli matrices and / is the 2x2 identity ma- 
trix. The components of the three-dimensional unit vec- 
tor n(k) are real functions of the wavenumber. 

Due to time reversal invariance of the Hamil- 
tonian e(— k) = e(k), £o(— k) = £o(k) and 
{n x (-k),n y (-k),n z (-k)} = {n x (k), -n y (k), n z (k)}. 
The energy spectrum of the two bands is £'±(k) = 
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£o(k) ± e(k) where the + and — signs correspond to the 
conductance and valance band, respectively. Since ft is a 
three-dimensional unit vector, it may be described by its 
azimuthal and polar angles which are also real functions 
of the wavenumber. 

n x (k) = sin i9(k) cos (p(k) 
^(k) = sin^(k) sin i^(k) 

n z (k) = cos^(k) 

Because of the time reversal properties of the Hamitonian 
i9(k) is an even function while </j(k) is an odd function of 
the wavenumber. The eigenvectors can be written as 

i?(k) 

|k ' +)= s -J?keMV 



The off-diagonal elements of the optical canductivity do 
not necessarily vanish and describe Faraday rotation if 
the symmetry group of the system is low enough. 

Applying (|10[) to chirally stacked multilayer graphene 
([T]) and taking into account both valleys, we obtain the 
well-known universal result 

M p 2 

muitiLG/, ,\ x n -n 

which is independent of the frequency-^. We note that 
the low energy model of multilayer graphene ([1]) is time 
reversal invariant only if both valleys are taken into ac- 
count. 

In the case of monolayer graphene the result (fTTj) was 
argued to be a consequence of the structure of the Dirac 
cones and the linear spectrum^. In the next section we 
show that the universal frequency dependence of the opti- 
cal conductivity can be found in a wide class of two-band 
systems. 



|k, 
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and are evidently determined by n(k). It is worth men- 
tioning that (k,±|<r|k,±) = ±n(k), i.e., the vector n(k) 
determines the pseudospin of the electrons and holes. 

At half filling photons excite electrons from the valance 
band to the conduction band. These transitions deter- 
mine the optical conductivity which can be obtained as 
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for positive frequencies. For negative frequencies 
Rea a p(uj < 0) = ReCT Q( g(|w|). Eq. © is valid only 
at zero temperature. For finite temperature one can give 
an analytical result only if £o(k) = 0. In this case ([5]) has 
to be multiplied by th(fej / 'AksT). Note that £o(k) does 
not appear in zero temperature conductivity. 

In tight binding approximation the Berry connection 
matrix element can be obtained as 



A_ + (k) =e 4v(k) 
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Substituting © into ([5]) and taking advantage of the 
time reversal properties of e(k), $(k) and ip(k) we found 
the following short expression for the optical conductiv- 
ity. For comparison, seei£. 
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This result shows that the optical conductivity of two- 
band systems can be calculated by identifying the vector 
n(k) and the dispersion e(k) in the Hamiltonian and then 
evaluating the sum over the Brillouin zone given in (1101) . 



IV. SCALE INVARIANT SYSTEMS 



The optical conductivity (fTU|) cannot be calculated 
generally. However, in this section we define a wide class 
of two-band systems for which the frequency dependence 
of the optical conductivity can be evaluated analytically. 

We use spherical coordinates in momentum space, i.e., 
k as the magnitude and {7} as the set of angle variables of 
the wavenumbcrs. If the dimension of the system is d the 
set {7} consists of d— 1 angle variables {71, 72, ■ • ■ , 7<j-i}- 
The gradient in momentum space can be written as 



dk 



d-l 



where efc and are unit vectors corresponding to the 
spherical coordinates k and 7^, respectively. Here Cj are 
coefficients of the Jacobian. Note that if the gradient 
acts on a function which depends on angle variables only 
the derivative is propotional to k^ 1 . 

In this section we focus on two band systems which can 
be described by the Hamilton operator with the following 
properties. We do not consider £o(k) because, as we have 
seen in the previous section, it does not influence the 
zero temperature optical conductivity. Additionally, we 
assume that the Hamiltonian depends on the magnitude 
of the momentum only through the dispersion and this 
dependence is power-law, i. e., e(k) = C({~f})k z where 
z is the dynamic exponent and C({j}) > 0. In this case 
n(k) = ^({7}) is a function of angle variables only. 

This assumption has lots of consequences. First, it 
follows that the conductance and valance bands touch 
each other at the Fermi level, i.e., there is no energy gap 
between them. Second, the system is invariant under di- 
latations which act on spacetime as r — > br and t — >• b z t 
where b > is the scaling parameter. This also means 
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that there is no characteristic energy or length scale in 
these systems. Third, the eigenvectors depend on angle 
variables only so the direction of the pseudospin is de- 
termined by the direction of the momentum only. Since 
this feature is a kind of generalization of the chirality of 
neutrinos and electrons in graphene, from now on, we 
may also refer to systems with the above characteristics 
as general chiral systems. Note that monolayer graphene 
and chiral stacked multilayer graphene are general chiral 
systems. 

The density of states in general chiral systems is a 
power-law function of the energy. 



9(e) 



(12) 



where the notation d {7} stands for the integration with 
respect to all angle variables and the Jacobian determi- 
nant is also incorporated. In mean-field theory systems 
with power-law density of states exhibit interesting criti- 
cal behaviour in quantum phase transitions accompanied 
by gap-opening^. 

In order to calculate the optical conductivity one has 
to determine the derivatives dn/dk a . Since ^({7}) is a 
function of angle variables only this derivative is propor- 
tional to k^ 1 . In the thermodynamic limit the summa- 
tion in Eq. (|10[) becomes an integral over the momentum 
space. This integral can be factorized into an integral 
with respect to angle variables and an integral with re- 
spect to k. Similarly to calculations presented in Ref.— , 
the latter integral can be carried out analytically in the 
case of general chiral systems and also determines the 
frequency dependence of the optical conductivity. 
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where k is the coefficient of the density of states (|T2"j) and 
K n fi is a dimensionless coefficient defined as 
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We note that in certain cases (for example in graphene) 
Eq. (|T3|) has to be multiplied by valley degeneracy. Since 
neither n nor K a p depend on the frequency, the opti- 
cal conductivity has universal uj( d ~ 2 ^ z frequency depen- 
dence. Using Kramers-Kronig relations the imaginary 
part of the conductivity can be obtained as Imcr(cj) ~ 
UJ ( d - 2 )/ z a t low frequencies. In a two-dimensional system 
(1131) provides a frequency- independent result. 



In one dimension there arc no angle variables so the 
unit vector n does not depend on the wavenumber which 
means that its derivative vanishes. It follows that in one- 
dimensional general chiral systems the optical conductiv- 
ity is identically zero for finite frequency. 

The result (fT3")) is universal in the sense that the mi- 
croscopic details do not influence the frequency depen- 
dence but are incorporated in prefactors only. The uni- 
versal behaviour can be understood on the basis of di- 
latation invariance. The optical conductivity scales as 
a(b,uj) = b 2 ~ d <j(l,b z uj) under dilatations defined above 
and the dilatation invariance demands da/db = 0. This 
equation results in cr(u>) ~ o/ d ~ 2 )/ 2 but says nothing 
about the prefactors. 

Dynamic scaling analysis is a very powerful tool to 
explore the frequency dependence of the optical con- 
ductivity in the vicinity of the critical point of phase 
transitions^—. Exactly at T c , power-law dependence 
with the exponent (d — z — 2)/z has been found^i. In 
quantum phase transitions one may use the effective di- 
mension d e ff = d+z instead of d which leads to the same 
frequency dependence as in Eq. (TT3"|) . It follows that 
in two-dimensional quantum critical systems the optical 
conductivity is also frequency-independent^. 

The scaling argument can be applied to other physical 
quantities. For example, the static polarization function 
of general chiral systems depends on the wavenumber as 
x(q) ~ q d ~ z in agreement with^i. Note that the scale 
invariance is closely related to the fact that there is no 
characteristic energy in the system. If an energy scale 
appears, for example due to gap-opening, Eq. (fT3| is 
not valid any more and the frequency dependence will be 
influenced by microscopic details, for example the way 
how the gap has been formed. 



V. CONCLUSION 

We presented a calculation of low frequency optical 
conductivity based on Kubo formalism. As a result, 
we obtained universal frequency dependence in systems 
which are invariant under dilatations and, hence, do 
not have characteristic energy scale. Scale invariance is 
closely related to the fact that the pseudospin is deter- 
mined by the direction of the momentum only. One may 
think of this property as some kind of chirality in a very 
general sense. The exponent of the frequency dependence 
is determined solely by the dynamic exponent and the 
dimension of the system. In two dimensions the opti- 
cal conductivity is independent of the frequency for any 
dynamic exponent in any scale invariant system. 
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Appendix A: Berry connection in tight-binding 
approximation 

The elements of Berry connection matrix arc defined in 
([5]). However, there is no analytical method to calculate 
uik, one has to apply some kind of approximation. 



In this appendix we investigate how the Berry connec- 
tion matrix can be expressed within tight-binding (TB) 
approximation. The localization properties of Wannicr 
functions has been studied extensively^. 

We restrict the Hilbcrt space to L atomic orbitals fb{^) 
(b = 1, . . . , L) per unit cell. Then one can define the 
Bloch functions 



ikR 



ip b {r-R) 



(Al) 



where R runs over all lattice vectors. We call this set 
of functions "natural" TB basis. Note that TB ba- 
sis is not necessarily orthonormal but can be Lowdin- 
orthonormalized^. We now assume that 

J dr^ 6k (r)V6'k'(r) = <5 kk/ <W . 

The Hamiltonian H(r) — — ft 2 A/(2m) + U(r) acts on 
as 

L 

ff(r) m (r) = fl«(kWk(r) + (H - PHP)<p bk (r) 

b' = l 

where H b , b (k) = J drip b/k (r)* H(r)ip bk (r) and P = 
X)b lyfekX^bkl is a projector onto the subspace spanned 
by the TB basis. In tight-binding approximation (H — 
PHP)ipb k (r) is neglected. 

In practice, the starting point in TB calculations 
is to determine H b / b (k). Then, the solutions of 
the Schrodingcr equation can be given as ifi k (v) = 
c/b(k)(/3fck(r) where the coefficients c;&(k) are deter- 
mined by the eigenvectors of H b i b {k). 

J2H b < b (k)ci b ,(k) =E lk ci b (k) 

b' 

The periodic part of fi k (r) can be written as 

u lk {r) = v / ?7e" lk Vzk(r) . 
Substituting in ([5]) we obtain 

A„,(k) = i^c i6 (k)V k c r6 (k)+ 
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w 



c /b (k)* Cw (k)^ e -' ikR 

R 



J drtp b {r - K)*rtpb'(r) . 



(A2) 

By introducing the notation |k, I) for the L-dimcnsional 
vector built up from the coefficients Qb(k) the first term 
12[) can be rewritten as i(k, /|V k |k, /). In the second 



term the matrix elements of the dipole operator appear 
which are taken between Wannier functions. Here we 
argue that these matrix elements are usually negligible 
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due to symmetry of atomic orbits or their well-localized 
behaviour but there are a few cases when it is not true. 
For instance, if both s and p orbits of an atom are taken 
into account in the tight-binding approximation the ma- 
trix element of the dipole operator is not negligible be- 



tween them. Nevertheless, in most cases (for example in 
graphene) the second term in (|A2|) can be neglected and, 
hence, 

A„,(k) = - / dTOjk(r)*Vk«i'k(r) w i(k, Z|V k |k, I) . 



